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Abstract 

The hemispherical 10" photomultiplier tube (PMT) R7081 from Hamamatsu Photonics K.K. (HPK) is 
used in various experiments in particle and astroparticle physics. We describe the test and calibration 
of 474 PMTs for the reactor antineutrino experiment Double Chooz. The unique test setup at Max- 
Planck-Institut fin* Kernphysik Heidelberg (MPIK) allows one to calibrate 30 PMTs simultaneously and 
to characterize the single photoelectron response, transit time spread, linear behaviour and saturation 
effects, photon detection efficiency and high voltage calibration. 
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1 Introduction 

Double Chooz [lj is a reactor antineutrino exper- 
iment located next to the nuclear power plant at 
Chooz, France. Its main physics goal is to mea- 
sure a non- vanishing value of the neutrino mixing 
angle #13 or to improve the current upper limit 
given by a global analysis [2\ 

sin 2 (2<9i 3 ) < 0.12 at 90% CL . (1) 

The two detector concept will substantially de- 
crease systematic errors and leads to a sensitivity 
of 

sin 2 (2<9i 3 ) < 0.03 at 90% CL (2) 

after 4 years of data taking [lj. Each of the 
two inner detectors of Double Chooz will be ob- 
served by 390 photomultiplier tubes (PMTs) to 
record light pulses produced by neutrino-induced 
reactions inside the gadolinium-doped scintillator 
target volume. The tasks of testing and cali- 
brating these PMTs were shared between German 
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Figure 1: The PMT Hamamatsu R7081 and the base circuit used for Double Chooz (taken from 0] and [3])- Rl, R2 = 
2.5 Mft; R3 = 180 kfi; R4 = 1.02 MQ] R5 = 1.5 MQ] R6 = 1 MQ] R7 = 499 kH; R8 = 300 kfi; R9 = 360 kQ; R10 = 430 
kfi; Rll = 680 kfi; R12 = 910 kfi; R13 = 750 kfi; R20 = 10 kfi; R21-R23 = 100 Q; R24 = 49.9 Q; C1-C3, C6 = lOnF; C4, 
C5 = 4.7 nF. 



General parameters of R7081 


Description 


Spectral Response 
Peak Wavelength 
Material of photocathode 
Effective photocathode area 
Structure of dynodes 
Number of dynodes 
Storage Temperature 


300 to 600 nm 
420 nm 
Bialkali 

220 mm dia. 

Box and Line 
10 

to +40 (°C) 



Table 1: General parameters of the photomultiplier tubes 
Hamamatsu R7081 [4]. 



groups (MPIK Heidelberg and RWTH Aachen) 
and Japanese groups (see Acknowledgements). 
Both groups worked with independent calibration 
setups. However, careful cross-checks between 
both setups showed all results in good agreement. 
This article describes the results obtained by the 
German groups, and the test setup which was built 
at MPIK Heidelberg. The Japanese setup and re- 
sults are described in |3J. 

The PMTs for the inner Double Chooz detector 
are 10 inch in diameter and have a hemispher- 
ical surface with a bialkali photocathode (type 
R7081MOD-ASSY) from HPK [4J (see Figure [l] 
and Table [TJ. They are made of glass which is op- 
timized for radiopurity. The specifications of the 
PMTs (see Table ^ have been validated for the 
purpose of a preselection for the Double Chooz ex- 
periment. Furthermore, detailed characterizations 
of the PMT behaviour are required for a precise 
understanding of the final detector response. 



These aims lead to the following calibration 
procedure which is performed for every PMT: 

1. investigation of the dark count rate during and 
after the stabilization phase 

2. determination of the high voltage (HV) value 
for a gain of 10 7 per photoelectron 

3. characterization of the single photon electron 
(SPE) response regarding spectral energy res- 
olution and peak to valley ratio 

4. exploration of relative transit times and the 
transit time spread for single events 

5. verification of the linear behaviour of the PMT 
signals up to high light intensities 

6. determination of the photon detection effi- 
ciency (product of quantum efficiency and col- 
lection efficiency) 

7. exploration of the afterpulse probability and 
their time characteristics 

8. search for flashing PMTs 

Besides the description of the test facility the fol- 
lowing sections focus on the analysis and the re- 
sults for the topics 1-6. The analysis of 7. and 8. 
is still ongoing and will be discussed in a separate 
publication. 

2 Test Facility 

The test facility for the above characterizations 
of the PMTs was built at MPIK in Heidelberg. 
The facility provides the possibility to measure 30 
PMTs simultaneously during one calibration-run. 
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Parameter 


Min. 


Typ. 


Max. 


Unit 


Supply voltage for gain of 1 x 10 7 


1150 


1500 


1650 


V 


Anode dark count 




4000 


8000 


1/s 


Peak to valley ratio at SPE 


2.5 


2.8 






Transit time spread (FWHM) 




3.4 


4.4 


ns 


After pulse probability (100 ns to 16 [is after main pulse) 






10 


% 


Quantum efficiency at 420nm 


22 


25 




% 



Table 2: PMT characteristics (and acceptance limits) at 25° C and supply voltage for a gain of 1 x 10 7 provided by 
Hamamatsu. 



During the whole calibration phase, three so called 
"reference PMTs" stayed in the setup and allowed 
to monitor the stability of the setup. Furthermore, 
the quantum efficiencies of these PMTs were ab- 
solutely calibrated by HPK which enables relative 
sensitivity measurements. In this section we de- 
scribe the setup, in particular the different light 
sources used, the electronics and the data acquisi- 
tion. 

2.1 Setup of the Test Facility 

The test facility was built inside a large 30 m 3 steel 
box. This, in the following called "Faraday room", 
is well protected against external electromagnetic 
fields and is closed light-tight. To prevent unde- 
sirable light reflections, the inner walls have been 
covered with black fire-proof rag. 

Figure [2] shows the topview of the Faraday 
room. The PMTs are mounted in slots of a racking 
system which consists of six rows and five columns 
(see Figure [3]). As denoted in Figure |2j each slot is 




LED fibre optics high voltage supply and 
signal line 



oriented towards the bis-MSB-ball, a central light 
source which allows one to illuminate all PMTs 
simultaneously. A metal arm at each slot feeds 
optical fibres coming from external light sources 



(see section 2.2) to illuminate each PMT individ- 
ually. The distance between the photocathode of 
a PMT mounted on the rack and the fibre attach- 
ment was chosen by two criterias with respect to 
the solid angle of the light spread at the end of 
the fibre. On one hand the photocathode should 
be fully illuminated, on the other hand light emit- 




Figure 2: Topview of the Faraday dark room containing 
the racking system for 30 PMTs. 



Figure 3: The PMT racking system with mounted PMTs. 
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ted from the fibre should not illuminate adjacent 
PMTs. Tests showed that both requirements were 
realized by a distance of about 50 cm between the 
end of the optical fibre and the photocathode. 

For further magnetic shielding, a /i-metal cylin- 
der covered by black paperboard is mounted in 
each slot. These cylinders have been produced for 
the use in the Double Chooz Experiment |5J. 

2.2 Light sources 

Three different light sources were used. 

• LED-boards: These VME controllable 12- 
channel boards were developed by the elec- 
tronic workshop at MPIK. The LEDs produce 
light at a wavelength of 380 nm with a light in- 
tensity that can be programmed for each chan- 
nel individually from the SPE range up to a 
few hundred PEs. The light is transmitted by 
optical fibres to each PMT. 

• Laser: For measurements which require a high 
precision in timing we used a Picosecond In- 
jection Laser (Advanced Laser Diode Systems 
A.L.S. GmbH) at a wavelength of 438 nm with 
a time jitter of 35 ps. Inside the light tight laser 
box (see Figure [2| the light pulses are widened 
with a diffusor lens and shines on a perforated 
plate where optical fibres are connected which 
lead to the PMT rack. Because of this config- 
uration, small differences of the light intensity 
are possible for different rack slots. 

• Bis-MSB-ball: This is a central light source 
of the test system. The glassblower work- 
shop at MPIK built a quartz ball with a di- 
ameter of 1 1 cm which was filled with a bis- 
MSB/ethanol solution. Bis-MSB is a wave- 
length shifter which is also used for the scintil- 
lator of the Double Chooz experiment [lj. The 
maximal spectral emission is at a wavelength of 
420 nm and matches the optimal sensitivity of 
the PMTs |4J. Four optical fibres are attached 
to guide LED signals via a quartz cylinder into 
the centre of the ball where the bis-MSB gen- 
erates an isotropic light emission to illuminate 
the full PMT array. 




Figure 4: The bis-MSB ball with attached optical fibres. 

2.3 Electronics and Data Acquisition System 
(DAQ) 

The front end electronics and DAQ system is lo- 
cated outside of the Faraday room. The PMT 
cables are fed through the Faraday room's wall 
and are connected to splitter boxes where the 
signals are decoupled from the supply voltage. 
The splitter boxes are designed and produced 
by the CIEMAT group (Madrid) for the Double 
Chooz experiment. The high voltage supply is 
a SY2527 Universal Multichannel Power Supply 
System (CAEN) j§] which provides a 500 mV res- 
olution for voltage setting and monitoring. The 
modules for processing and acquiring the PMT 
signals are based on standard NIM modules and 
a VME system, respectively. Figure [5] shows a 
schematic illustration of the data acquisition sys- 
tem. 

In order to exploit the full dynamic range of 
the used acquisition modules and to provide the 
required logical signals, several NIM- modules are 
used: 

• 4x sixteen channel amplifier (fixed gain xlO, 
Phillips Scientific (PS) [7J Mod. 776). 

• 4x octal discriminator (threshold > 10 mV, PS 
Mod. 710). 

• 4x octal linear Fan-Out (PS Mod. 748). 
The used acquisition modules are listed below: 

• 32 channel charge-to-digital-converter (QDC 
V792, 12 bit resolution (400 pC), CAEN [§]) 
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• xl5 




• xl5 



Splitter | 



Splitter | 



Splitter 



Splitter 



HV supply 



Laser 
LED boards 



Ampl. xlQ |-| Ampl. xlp] 



Ampl. xlQ ^ Ampl. xlO 



Ampl. xlQ ^ Ampl. xlp] 



Ampl. xlQ |-| Ampl. xlp] 
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TDC 
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III 


TB 1 


TB 2 



MPIC board 



Figure 5: Schematic diagram of the electronics setup and the DAQ. 



• 32 channel time-to-digital-converter (TDC 
V775, 12 bit resolution (400-1200 ns), CAEN) 

• 2 x 16 channel Scaler (V260E, 24 bit channel 
depth / 100 MHz counting rate, CAEN) 

• The Double Chooz Trigger Board (DC-TB), 
RWTH Aachen [8]. 

All trigger signals, gate widths etc. for the mea- 
surements were provided by the MPIC-board, 
which was developed at the electronic workshop at 
MPIK. This board provides up to 15 logic signals 
with free programmable gate widths and delays. 

3 Measurements and results 

3.1 Dark count rate behaviour 

Dark count events are caused by thermal emis- 
sion of electrons from the photocathode which are 
not distinguishable from events caused by regular 
photoelectrons. The major part of these events 
are caused by single electrons. We measured the 
dark count rate with a discriminator threshold at 
1/4 SPE pulse amplitude height (corresponding to 
1 mV for an unamplified signal at a gain of 10 T ). 

If PMTs are exposed to light, which is unavoid- 
able during the installation of the PMTs at the 
test setup, the dark count rate is increased during 



the first hours after turning on the high voltage. 
This phase of stabilization was monitored for a 
period of typically 20 hours with a scaler and the 
DC-Trigger-Board (see Figure |6| left). After this 
period of time, the rate has decreased to a reason- 
ably stable value and all calibrations which will 
be described in the following sections were per- 
formed. Afterwards, the dark count rate was ob- 
served again for about one day to determine the 
final value and to inspect its stability. 

The dark count rate depends on the temper- 
ature. During the measurements the tempera- 
ture inside the laboratory at MPIK was stabi- 
lized by air conditioning at (20 ± 1)°C. We ob- 
served a change of the rate of roughly 100 c/s 
(counts/second) per 1°C temperature difference. 
Inside the Chooz underground laboratory the tem- 
perature is about 14° C, therefore the dark count 
rate is expected to be lower in the final detector. 

The specification for the Double Chooz experi- 
ment is a dark count rate below 8000 c/s per PMT. 
With one exception, all tested PMTs showed much 
smaller values with an average of 2200 c/s at 20° C 
with a standard deviation of 500 c/s. 

To quantify the stability of the dark count rate, 
we calculated the root mean square (RMS) of 
taken data points during the stable phase (see Fig- 
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5000 




Figure 6: Left: typical dark count rate behaviour during the first hours after installation. Right: example for an unstable 
dark count rate. 




200 300 400 

arbitrary photomultiplier number 



Figure 7: RMS value from the last 20 minutes in the dark 
count rate measurement. 



ure[6j right) to compare it with the value which is 
expected for statistical fluctuations. If the RMS 
was significantly higher (see Figure [7]) the PMT 
was not used for the Double Chooz detector. 

3.2 Gain Calibration 

For this measurement, we used the LEDs at SPE 
level. The light intesities were adjusted in such 
a way that on average only one out of ten trig- 
gers gives a PMT signal. Assuming that the num- 
ber of created photoelectrons per event npE is 
poisson distributed, the probability of npE > 2 
is only about 0.5%. Therefore almost pure SPE 
spectra were observed. The total trigger rate of 
LED pulses was 1 kHz. For each PMT, we took 
nine charge spectra (over 400 s each which cor- 



responds to 400k events) at high voltage values 
HVj (j = 1, . . . , 9) ranging between ± 200 V 
around the nominal high voltage HV nom (which is 
the value for a gain of 10 7 as provided by HPK). 
The signal current arriving at the QDC was inte- 
grated within a time window of 200 ns. This leads 
to charge spectra for which an example is shown 
in Figure [8] To analyse these spectra the pedestal 
peak was fitted by a Gaussian Ped(i) 



Ped(i) = 



Ped 



27TCrped 



exp 



{i ~ MPed)' 



(3) 



To describe the SPE peak a combination of an 
exponential and a Gaussian SPE(i) was used 



SPE(i) = / ^ SPE exp 



27RT S PE 



(i - ^spe) s 
2ct 2 



iV] 



Exp 



exp 



(4) 



7Vp ec i and A/spe are the number of entries in the 
pedestal and the SPE peak, respectively. /ip e d and 
Mspe are the mean values of the fitted Gaussians, 
<jspe and aped their standard deviations, A^Exp 
and r are the free fit parameters of the exponen- 
tial which is an empirical function to describe bad 
amplified events, i is the number of the QDC-bin 
and corresponds to a charge of i- 400/4096 pC, i m i n 
is the endpoint of the pedestal fit and the start- 
ing point for the exponential fit which is roughly 
determined by the valley position. 
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Figure 8: SPE charge spectrum in linear and logarithmic view, together with fit functions for pedestal and SPE peak. 
The red line shows raw QDC data. The pedestal peak was fitted by a Gaussian (green). The functions whose combination 
describe the SPE peak (see text) are shown in blue and orange. The black line is the combination of all fit functions which 
describes the complete spectrum. 



The fit functions and their combination are 
drawn in Figure [8] showing model and data in good 
agreement. From the results of both fits the av- 
erage charge Q corresponding to a SPE-signal is 
calculated 

Q = MSPE - MPed • (5) 

For the nine different high voltage values HV J5 
the corresponding charges Q(HVj) are shown in 
Figure [9j The Q(HV) dependence is fit with a 
power law function 



Q(HV) = 8.0pC- 



HV 
HV ( 



opt 



(6) 



The free fit parameters are HV opt and a, where 
HV op t is the searched value for a gain of 10 7 which 
corresponds to an average charge of 8.0 pCQ 
The obtained results of HV op t for all tested PMT 
are presented in a histogram shown in Figure [To| 
(left), together with the values HV nom provided by 
HPK. The good correlation between these values 



is shown in Figure 10 (right). 



For the verification of the obtained results and 
a detailed characterization of the SPE response 
of each PMT an additional SPE charge spectrum 



1 A gain of 10 7 corresponds to a charge of 10 7 e = 
10 7 • 1.602 • lCr 19 C = 1.602 pC. An additional 50ft back 
termination of the PMTs causes a division by a factor of 2, 
finally an amplification ( x 10) in the electronics chain leads 
to a value of 8.0 pC. 




1200 1250 1300 1350 1400 1450 1500 1550 1600 1650 

HV in V 

Figure 9: Example of a power law fit for determination 
of HVopt- The black points represent the charge Q which 
was determined at a certain high voltage value according to 
equation |5j. The solid line is the fitted power law function 
equation J6l. 



was recorded at the previously determined HV op t 
value. The charge value for this spectrum is deter- 
mined with the above procedure and is compared 
to the expected value of 8.0 pC. The results for 
all PMTs are shown in Figure [TT] The gain values 
deviate from the expected gain with a Gaussian 
distribution and relative RMS of 1.3%. 

Finally, we investigate the stability of the three 
reference PMTs which are repeatedly measured in 
each data run at fixed positions in the rack. Figure 
12 shows the determined HV opt values as a func- 
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Figure 10: Left: Distribution of determined values HV op t (solid line) for a gain of 10 7 compared to the values HV n0 m 
provided by HPK (dotted line) of all tested PMTs. Right: HV op t values versus the corresponding values HV nom - Both 
design and actual HV values are confirmed to be consistent. 
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Figure 11: Average charge value of SPE signals at HV, 
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Figure 12: HV op t of the reference PMTs versus the run 
number. 



tion of the consecutive run numbers of successive 
PMT calibrations (corresponding to a period of 
about three months in total). The obtained val- 
ues HVopt agree within a few Volts, however a 
peculiar step is observable at run 7. Before this 
run, it was necessary to change some amplifier and 
fan-out channels which required re-calibrations for 
these channels. The visible step is consistent with 
an assumed systematic error of the calibrations of 
about 1% (which was determined using a pulse 
generator). In Figure 13 the determined pedestal 
values are shown in a histogram. The average rel- 
ative spread of less than a half percent of the de- 
termined pedestal values shows a good stability of 
the test setup. 



3.3 SPE response - charge and transit time 

The PMT response for SPE pulses was character- 
ized in detail using the high voltage HV op t- The 
DAQ provides the possibility to record simulta- 
neously charge and transit times for each event. 
A typical SPE charge spectrum is shown in Fig- 
ure [8] The charge resolution of the SPE peak is 
typically AE w 2pC, or to give a relative num- 
ber AE/E « 0.25. The peak-to- valley-ratios P/V 
(defined as the ratio between the maximum of the 
SPE peak and the minimum in the valley between 
the pedestal and the SPE peak) is in the range of 
3.2 < P/V < 5.5. The results for all tested PMTs 



are shown in Figure 14 These numbers were im- 
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Figure 13: Histogram of the determined pedestal values 
of the reference PMTs. 
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Figure 15: Results for all tested PMTs regarding transit 
time spread. 




PMT Serial No. 

Figure 14: Results for all tested PMTs regarding Peak- 
to- valley-ratio. 



tribution: pre- and late-pulse^] Pre-pulses are 
appearing about 30 ns before the main peak. They 
are presumably caused by photons which pass the 
photocathode and hit the first dynode producing a 
photoelectron. In this case one step of the electron 
multiplication is missed, and the charge of pre- 
pulses is lower than the SPE value. Late-pulses 
appear if photoelectrons are (in)elastically scat- 
tered at the dynode, reflected and re-accelerated. 
These late-pulses appear up to 80 ns later than 
the main pulses and have a slightly smaller charge 
than regularly accelerated photoelectrons. The 
charge-transit time correlation is shown in Figure 
16 (right). The probability of pre-and late-pulses 
are roughly 0.1% and 3%, respectively. 



proved substantially from previous smaller values 
after the introduction of the /i-metal shields. 

The absolute transit times were measured in- 
cluding an additional offset related to the electron- 
ics and the used light source. The distribution of 
relative transit times is shown in Figure [l6] (left). 
Most of the events occur at roughly the same time 
after the time of the Laser trigger which defines 
t = 0. This main peak can be described by the 
sum of two Gaussians. The transit time spread 
is defined as the width of this peak (full width 
at half maximum), the results of all tested PMTs 



are shown in Figure 15 The average transit time 
spread (FWHM) is 2.8 ns with an RMS of 0.2 ns. 

In addition to the main peak two more classes 
of events can be identified in the transit time dis- 



3.4 Verification of linear behaviour up to high 
light intensities 

The expected light levels of neutrino events in 
the Double Chooz experiment are typically in the 
range of a few PE. Nevertheless, signals induced 
by cosmic muons, fast neutrons etc. may exceed 
the SPE region substantially. Therefore we stud- 
ied the linear behaviour of the PMTs up to light in- 
tensities corresponding to 30 photoelectrons (PE) 
(in a few cases, up to 300 PE). 

The measurement proceeds similar to the gain 
calibration except for the use of the bis-MSB-ball 
as central light source. Because of the higher 



2 Please note that the used single-hit TDC is not sensitve 
for afterpulses (caused by ionized residual gas) since only 
the first event after the light trigger is recorded. 
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Figure 16: Left: typical distribution of relative transit times including main-, pre- and late pulses. Right: charge versus 
transit time plotted event by event. 



light intensity almost every pulse of the bis-MSB- 
ball leads to PMT signals. Therefore, additional 
DAQ triggers without light pulses were taken for a 
precise pedestal determination. A typical charge 
spectrum is shown in Figure [T7[ The pedestal was 



fitted with a Gaussian (see equation (J3|), and the 
spectrum was corrected for this pedestal. The 
charge is converted into the corresponding num- 
ber of PEs, using the values determined during the 
gain calibration (equation Q). The expected dis- 
crete Poissonian distribution with the mean value 
ji is smeard by the SPE-resolution <j. Therefore 
a convolution of the Poissonian with a Gaussian 
was used to fit the multi-PE data (MPE) and to 
determine /i 




Number of PE 



Figure 17: Typical charge spectrum for light intensities 
in the multi-PE region. The data (red) were fitted by the 
function described in equation (JT]) (black). 



MPE(i) 



k=l V 



k\ 



exp(— fi) x 



exp 



(i - kf 
2ka 2 



(7) 

The free fit parameter is \i which is the average 
number of photoelectrons. N is the total num- 
ber of entries and <j the SPE resolution, i is the 
number of the bin of the QDC corresponding to a 
charge calibrated in PE. 

For each PMT, nine charge spectra for the same 
HVj values as for the gain calibration in section 
3.2 are measured. The determined fij are con- 
verted to the corresponding charges Qj which are 
plotted as function of the HVj values in Figure 



18] A power law function is fitted to the data and 
compared to the gain calibration (equation (|6|) 



Q(HV) = 7-8.0 pC- 



HV 



HV, 



opt 



(8) 



The free fit parameters are the relative number of 
PE for this light level I and the exponent a. HV op t 
had been determined earlier during the gain cal- 
ibration (see section 3.2). In Figure 18 both HV 



scans are shown in double logarithmic scale. The 
slope QfMPE for higher light intensities is similar 
to the slope aspE at small light levels indicating 
a good linear response. A good linear behaviour 
should lead to a ratio K = QMPE close to unity. 
The resulting values of K for all tested PMTs are 



shown in Figure 19 The mean of the distribution 
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Figure 18: Average charge as a function of correspond- 
ing HV values at two different light intensities in double- 
logarithmic scale. 
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Figure 19: Distribution of the K ratios of all tested PMTs 
and fitted Gaussian. 



is close to 1 with a spread of only 0.6%. No indi- 
cations of non-linear behaviour up to signals cor- 
responding to 30 PE were observed for any PMT. 

For a few PMTs, further studies of the linear re- 
sponse have been performed after the completion 
of the calibration measurements with even higher 
light intensities. In order to avoid saturation of 
the QDC one amplifier in the DAQ chain has been 
removed. The results lead to the conclusion of a 
good linear response up to 250-300 PE when op- 
erating the PMTs at HV opt . 

3.5 Determination of the photon detection effi- 
ciency 

As an important characteristic of the PMTs, we 
determined the photon detection efficiency e in 
order to identify bad PMTs. The efficiency e is 
defined as the product of the quantum- and the 
collection efficiency e = QE x CE. The quantum 
efficiency (QE) indicates the ratio of the number 
of created PEs and the number of photons hitting 
the photocathode. The collection efficiency (CE) 
determines the number of PEs which reach the 
first dynode. With our setup it was not possible to 
measure both quantities (QE, CE) independently, 
however from an experimental point of view only 
the product e is of interest. Following the infor- 
mations given by the manufacturer, the collection 
efficiency is 0.9 and the quantum efficiency about 
0.25 in the region of interest at a wavelength of 
420 nm. The photon detection efficiency s de- 



pends on the voltage between the photocathode 
and the first dynode, the structure of the electri- 
cal field as well as on the effective area of the first 
dynode. 

The bis-MSB-ball is used to illuminate the 
complete PMT array and the absolute number of 
photons reaching the cathodes is unknown. There- 
fore, the performed calibrations are relative, com- 
paring each PMT with the quantum efficiencies of 
three absolutely calibrated reference PMTs pro- 
vided by HPK. 

Because the light source is not perfectly stable 
in time, the light intensity of each event was mon- 
itored by a fixed reference PMT which remained 
in the same slot during all measurements. Fur- 
thermore, the geometric dependencies were cali- 
brated by moving two additional reference PMTs 
subsequently slot by slot throughout the complete 
array. For each slot s a geometry factor g s was 
determined, using the ratio of the mean values of 
the fixed PMT /i fix and the reference PMT at the 
corresponding slot /i s , defining g s := fin x /fi s . For 
both reference PMTs (ref = 1,2) the procedure 
was repeated for two different light intensities Ii^- 
The final calibration for each slot was obtained by 
averaging 

By comparing the response fipMT of a PMT at slot 
s with the simultaneous response /ifi x of the fixed 
PMT it is possible to evaluate the value £pmt of 
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Figure 20: Photon detection efficiency e = CE x QE of 
all tested PMTs and fitted Gaussian. 
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Figure 21: Photon detection efficiency e versus the corre- 
sponding ScB-value. 



the unknown PMT relative to the reference PMTs 
used during the geometric calibration 



£PMT 



MPMT re f 

y s 

Mfix 



£ref • 



(10) 



To increase statistics, the measurements were 
done at three different light intensities in the 
multi-PE region. In addition (as mentioned 
above), we used two independent reference PMTs 
to reduce systematic errors regarding the refer- 
ence QE and in order to have a backup solution if 
one of the reference PMTs should show any unex- 
pected irregularities (which did not happen). As 
final value we used the arithmetic mean of the six 
results (three intensities and two reference PMTs) . 
For all charge spectra the same analysis methods 
were used as described in section [3741 

The geometrical calibration procedure which 
leads to equation ([9| performed at the begin and 
was repeated twice during the calibration phase to 
control the stability of the system. It turned out 
that the used light source (the bis-MSB-ball) was 
not perfectly stable with respect to geometrical 
anisotropics. This means the values g v s ei were af- 
fected by changes in time, either slowly or stepwise 
which could not be determined by a backdated 
analysis. These shifts were probably caused by 
small changes of the position of the optical fibers 
leading light from the LED boards inside the crys- 
tal ball. Analysing the response of the two refer- 
ence PMTs with respect to the fixed PMT (all 
these PMTs were present during the complete pe- 
riod of calibration) and comparing these values 



with the results of the three geometrical calibra- 
tions we estimated the relative systematical error 
of the photon detection efficiency to As < 5%. 
This precision is lower than expected before the 
calibration campaign, however it is sufficient to 
identify bad PMTs which was the main goal of 
these measurements. In total, there were only 
3 PMTs with a photon detection efficiency lower 
than 19% which were not dedicated for installa- 
tion in the Double Chooz detector. 

HPK provided a value describing the sensitiv- 
ity for blue light (ScB = sensitivity cathode blue 
index) for each PMT, which is expected to be cor- 
related to the determined efficiency values e. In 
Figure [21] the results e of all PMTs are plotted 
versus their ScB-value. The correlation is obvi- 
ous, although not very strong which is caused by 
the uncertainties of both values and differences in 
the collection efficiencies which are included in e. 



4 Conclusions 

For the Double Chooz experiment, 474 photomul- 
tiplier tubes were tested and calibrated regarding 
the dark count rate, gain as a function of high 
voltage, SPE response, transit time spread, linear 
behaviour and photon detection efficiency. One 
PMT showed no signal at all, another one lost its 
vacuum after the calibration measurements for an 
unknown reason. For 5 PMTs we noticed an un- 
stable dark count rate. In addition, one of them 
showed an inacceptable dark rate of more than 
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10000 c/s . For 3 PMTs we determined an pho- 
ton detection efficiency lower than 0.2. In total 8 
PMTs are not dedicated to be installed inside the 
Double detectors and will be kept as spare and to 
perform additional offsite tests in the future. 
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